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We show how to apply the macroscopic fluctuation theory (MFT) of Bertini, De Sole, Gabrielli, Jona-Lasinio, 
and Landim to study the current fluctuations of diffusive systems with a step initial condition. We argue that 
one has to distinguish between two ways of averaging (the annealed and the quenched cases) depending on 
whether we let the initial condition fluctuate or not. Although the initial condition is not a steady state, the 
distribution of the current satisfies a symmetry very reminiscent of the fluctuation theorem. We show how the 
equations of the MFT can be solved in the case of non-interacting particles. The symmetry of these equations 
can be used to deduce the distribution of the current for several other models, from its knowledge [29] for the 
symmetric simple exclusion process. In the range where the integrated current Qt ~ \ft, we show that the 
non-Gaussian decay exp[— Q\ /t] of the distribution of Qt is generic. 
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I. INTRODUCTION 

The study of the fluctuations of currents of energy or of particles is central in the theory of non-equilibrium systems. Over the 
last decade, the macrosopic fluctuation theory (MFT), a theory of diffusive systems maintained in a non-equilibrium steady state 
by contact with two heat baths or two reservoirs of particles, has been developed [1, 2]. This theory was first implemented to 
give a framework to calculate the large deviation functional of density profiles in non-equilibrium steady states 10, SlaS la. Si- 
It was then understood that it could also be used to predict the distribution of the current through non-equilibrium diffusive 
systems 

The macroscopic fluctuation theory gives a large scale description of lattice models such as the symmetric simple exclusion 
process (SSEP) or the Kipnis Marchioro Presutti model fl6l [l7l [l8i [l9Tl . At the microcopic level, several properties of these 
models can be obtained using numeric al llT/l[l8l,l2"oll . perturbative fell 12211 . or exact approaches ll23ll . such as the matrix method 
ll24l EM l26ll or the Bethe ansatz lfl5i l27l l28ll . Whenever the comparison has been possible, it is remarkable that a perfect 
agreement has been found between the results (on the large deviations of the density profile luil l25i l26ll or on the probability 
distribution of the current 

COE US, IH) obtained by these microscopic approaches and the predictions of the macroscopic fluc- 
tuation theory (5l [3, [l0|] . Moreover the MFT led to the prediction of rather surprising properties of diffusive systems, such as 
the possibility of phase transitions lfl2l Il3l [Till in the large deviation function of the current, or the universality of the cumulants 
of the current on the ring geometry lll5ll . So far the MFT has only been used on systems at equilibrium, or in non-equilibrium 
steady states. 
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FIG. 1 : The step initial condition with a density p a at the left of the origin and pt at the right of the origin 



In a recent work [29], we considered the fluctuations of the integrated current Q t through the origin of the SSEP, starting with 
a non steady state initial condition : a step in the density profile at the origin with density p a on the negative axis and density 
Pb on the positive axis, as shown in figure Q] The SSEP is one of the simplest lattice gas models, and has been extensively 
studied in the theory of non-equilibrium systems lf30l l3ll l32ll . The distribution of the integrated current Qt, for the SSEP, is 
related to the time decay of constrained one dimensional Ising models B33I1 . In the SSEP, particles diffuse on the lattice with 
nearest neighbor jumps and a hard core interaction which enforces that there is never more than one particle on each site (in 
practice, the configuration at time t is specified by a binary variable n (t) — {1 or 0} on each lattice site, which indicates whether 
site i is occupied or empty; the dynamics is such that these occupation numbers are exchanged at rate 1 between every pair of 
neighboring sites on the lattice). Using the Bethe ansatz and several identities proved recently by Tracy and Widom 11341 l35ll3al 
for exclusion processes on the line, we were able to show 1I29I1 that the generating function of the total flux Q t of particles 
through the origin during a long time t takes the form 



with/i(A,/9 a ,P&) given by 



n(X,p a ,Pb) = - dklog 

7T ' 



1 + use" 



(1) 



(2) 



and where oj is a function of p a , pb and A 

lu = p a {e x - 1) + p b {e- x - 1) + p aPb {e x - l)( e - A - 1) . (3) 
Beyond the fact that that /i is a function of the single parameter u, which was proved in [29], one can see from (111213b that 

1 . All the cumulants of Q t grow like y/i. 

2. p(X, p a , pb) satisfies a symmetry very reminiscent of the fluctuation theorem ll40l [4lL [42l l43l [441. [45ll : 



p a ,Pb) = V ^-A + log- 



Pb 



Pb 



log 



Pa 



1 - Pa 



Pa,Pb 



(4) 



(this is because w in © is left unchanged by this symmetry). 



3. For technical reasons in the way that (O was derived in [29], we had to impose the condition that \u\ < \j2j~K. If one 
assumes that the range of validity of <(2J extends to all u > — 1, one gets that p ~ ^(logw) 3 / 2 for large uj, which would 
imply that for large q 



Probability 



(3H 



exp 



12 



Vt 



exp 



12 t 



31 



(5) 



The goal of the present work is to see how the above results <JT][5j, obtained for the SSEP with the step initial condition of 
figureQ] can be understood from the point of view of the MFT and how they can be extended to more general diffusive systems. 

When the dynamics is stochastic, the integrated current Q t through the origin depends both on the history (i.e. on all the 
updates between time and time t) and on the initial condition (which, for the SSEP, is drawn according to a Bernoulli measure 
of mean p a on the negative axis (i < 0) and pb on the positive axis (i > 1)). Very much like in the theory of disordered systems, 
where one can distinguish between an annealed average (where the partition function is averaged over all the realizations of the 
disorder) and a quenched average (where the partition function is calculated for a typical realization of the disorder), one can 
define here two expressions of p(\) : 
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• the annealed case where, as in the derivation of (Q][3]) in l29ll . one averages e A( ^* both on the history and on the initial 
condition 



Ma„n„ealcd(A) = Hm log [(e AQf ) history , initial 



I condition 



(6) 



It turns out that the initial conditions which dominate the average are atypical as shown in figure |2] 
• the quenched case, where one averages e A< 5* only on the history for a typical initial condition 

AVc„chcd(A) = lim \ ( log [<e AQt ) histo J )...,,. (7) 

t^oo y/t \ I History j / ; nitia i con dition 

The difference between these two averages, and their influence on the distribution of the current, has already been studied for 
the totally asymmetric exclusion process (TASEP) using the microscopic dynamics 13711 . 

In section HU we formulate the calculation of both /i a imcaicd and ^(quenched in the framework of the MFT. In section [TTH we 
see that /^annealed satisfies the symmetry for general diffusive systems, and for general non-steady state initial conditions. 
No such symmetry seems to hold in the quenched case. In section HVl we consider the case of non-interacting random walkers, 
where both jU ann eaied an d /^quenched can be determined exactly. In section |VJ we show that, for the SSEP, the single-parameter 
dependence (fJJ of /J, anne aied can be understood from a remarkable invariance of the MFT. In section I VII we obtain bounds on 
the decay of the distribution of Q t which shows that (|5]l is generic for a broader class of diffusive systems. 




FIG. 2: Average rescaled density p(x, r) (see J14b ) for the SSEP, in the annealed and quenched cases, when A = 1.5, p a = 0.6 and pt = 0.4. 
While the initial profile is a step function in the quenched case, it deviates from it in the annealed case. 



II. THE ANNEALED AND THE QUENCHED AVERAGES 

In this section we show how the macroscopic fluctuation theory U,^] can be used to calculate the generating function of the 
integrated current Q t when the initial condition is a step density profile. The theory is in principle valid for arbitrary diffusive 
systems with one conserved quantity, such as the number of particles or the energy. Here, for simplicity, we consider the case of 
a one dimensional lattice gas where a configuration is characterized by the numbers rii of particles on each site i. 
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Imagine first that this one dimensional system has a finite length L, and that it is in contact, at its two ends, with two reservoirs 
of particles at density p a and p . In this finite geometry, the system's stochastic evolution reaches a steady state, where the flux 
Qt of particles during a long t has a certain average (Q t ) and a certain variance (Q 2 ) c = (Qt) ~ (Qt) 2 - Close to equilibrium, 



i.e. when the densities of the two reservoirs are close (p a 
expects II 1 CH. 13211 that 



Pb 



r with p a — Pb r"), and for a large system size L, one 



lira 



(Qt) „ D(r) 



(Pa - Pb) 



(8) 



and 



lim 



(Q 2 t)c a(r) 



(9) 



where D(r) and er(r) are two functions which characterize the transport of particles through this diffusive system. 

At equilibrium (for p a = Pb = r ), the weights of all microscopic configurations are given by the Boltzmann weights. For 
large L, if one introduces a rescaled position < x = i /L < 1, the probability of observing a given density profile p(x), when 
the two reservoirs are at the same density r, satisfies l32ll 

Pro cq .(p(a;)) x exp[-LJ" oq .(/o(x))] , 



where the large deviation function JF oq . is given by 



l r 



f( p (x))-f(r)-(p(x)-r)f(r) 



dx . 



and f(r) is the free energy per site of the equilibrium system at density r (defined as f(r) - 
Z(N, L) the partition function of the system with N particles on L sites). One can show ]3C 
theorem, which is satisfied at equilibrium, implies that 



f"(r) 



2D(r) 
o-(r) 



(10) 

lirxii^ooClog Z(Lr, L))/L) far 
',] that the fluctuation dissipation 



(11) 



For a diffusive system on a one dimensional lattice of L sites, in contact with two reservoirs at densities p a and p^, the average 
density (iii(t)) near position i and time t, and the total flux of particles Qi(t) through position i between times and t, are 
expected to follow diffusive scaling laws. For large L, and for times of order L 2 , they take the form 



{rii{t)) = p[~,-L 



V L 1 



i t 

and Qi(t) ^ L( l[j j ^Jj2 



From the large deviation hydrodynamics theory lfl4ll30l[32ll38i[39ll . the probability of observing a certain density profile p (x, r) 
and current profile j (x,t) = dq(x, r) jdr over the rescaled time interval < t < t/L 2 is expressed as 



Pro 



(ip{x,T~),j(x,T)} 



exp 



rt/L 2 j.1 
-L dr dx 

Jo Jo 



j( x ,r)+D(p(x,r))^ 



2ct(p(x,t)) 



(12) 



where D(p) and a(p) are defined as in d8l9t . Expression ( fT2l simply means that, locally, Fick's law (j = —D(p)p') is satisfied 
everywhere up to Gaussian current fluctuations of variance cr(p). The conservation of the number of particles, 
rii(t) — rii(0) = Qi-\(t) — Qi(t), becomes a conservation law on the rescaled density and current profiles : 



dp 



dx 



(13) 



For a non-steady state initial condition, as in figureQ] the system size is infinite. If one observes the fluctuations of the current 
over a long time t, one can introduce a characteristic length \ft. The average density (rij(i')) near site i and the integrated 
current Qi(t') between times and t' < t then become scaling functions of the form 



M{t')) 



= n - 



t' 



and Qi(t') = Vtq 



(14) 
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and the probability of observing such rescaled density and current profiles is given by 



Pvo^{p(x,T),j(x,T)}^ xexp 



-Vi [ dr [ 

Jo J-t 



dx 



dp(x,r) 



j( x ,r) + D(p(x,r))^ 



2a(p(x,T)) 



(15) 



The integrated current Qt through the origin during time t can then be written as 



roc 

Qt=Y\ rn(t) - - V* / dx [p(x, 1) - p(x, 0)] . 



(16) 



Moreover, when the initial condition is a local equilibrium configuration at density p a on the negative axis and density pf, on 
the positive axis, as in figureQ] the probability PrOi n iti a i of the initial profile p(x, 0) is given by 



ProinitiaifXz, 0)) >c exp -Vi Fmit.(p(x, 0)) 



where (fToTTB 



^initXP^) )) 



f(p(x,0)) - f(r(x)) - (p(x,0) - r(x))f(r(x)) 



p(x,o) 2D(z) 
dx I dz (p(x,0) -z)- 

r(x) 



dx 



and r(x) is the step density profile 



(0(x) is the Heaviside function). 



a{z) 

r{x) = (1 - 0{x))p a + 9(x)p b 



(17) 

(18) 
(19) 

(20) 



The annealed case 

Therefore d6l 1 151 161 [TBI lead to the following expression for /i an nnealed 



Mannnealed(A) = max < -T ia it. (p{x, 0)) + A / dx [p(x, 1) - p(x, 0)] - / dr dx 

p(x,T),j(x,T) | Jq J Q 



dp(x,r) 



2a(p(x,r)) 



(21) 



Finding the optimal p(x, r) and j(x, t) in (|2TT > has to be done carefully because they are related by the conservation law ( [Oi l 
As shown in Appendix[A](see also [5, 8]), one can replace the variational form d2lT > by another variational form : 



Mannnealed(A) = max <^ -^init. (p(x, 0)) + A / dx [p(x, 1) - p(x,0)] 

p (x,t),h(x,t) y J Q 



dr / dx 



.dp{x,r) dH(x,r) dp(x,r) a(p(x,r)) ( dH(x,r) 
H ( x > T > fl _ + D (P( X > T >> HI 2 I 



<9x 9a; 



<9;C 



where p(a;, r) and H(x, r) are independent. The optimal p(x, r) and H(x, r) in (l22l ) satisfy 



dp(x, t) d 
dr dx 

dH(x, t x 



D(p(x,r)) 



dpjx, t) 
dx 



d_ 



(x(p(x,T)) 



8H(x,t) 



dx 



dr 



= -D(p( X) r)) 



2 H(x,t) <j'( P (x,t)) fdH(x,r) 



dx 2 



V dx 



with the boundary conditions 



#0,1) = A0(x), 

H(x,0) = X9(x)+2 / —rrdp, 

Jr(x) Cr(p) 



(22) 

(23) 
(24) 

(25) 
(26) 
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where we have used ( fTTT i and ([T8~l l. 

Thus in annealed case one can use either d2"TT i or d22l|26l l to obtain /i an nncaicd(A). Using the fact that p(x, t) satisfies d23l ) and 
that p{x, t) and H(x, t) have limiting values (p a , p ) and (0, A) as x — ► ±00, one can simplify (f22b to get 



Man„nc al cd(A) = -.F init . 0)) + A f°° dx[p(x , 1) ~ p(x , 0)} - f dr f°° dx ^^ ( ~ ) (27) 



The quenched case 

In the quenched case, the main difference is that p(x, 0) is no longer allowed to fluctuate. Therefore the boundary condition 
i at r = is replaced by 



p{x,0) = r(x) , 



and ( f2TT > becomes 



/"quenched (A) = max < A / dx \fi(x, 1) - p{x,0)] 



dr / dx 

— oc 



dp(x.r) 



2a(p(x, T )) 

with the max taken over all the profiles satisfying and d28l >. In terms of the field H, one gets 

poo pi 

/iquenchcd(A) = A / cfe[p(a;, 1) - p(x, 0)] - dr 
Jo Jo 

with p and iJ satisfying d23ll25l ) but with ([28]) instead of 



dx 



(28) 



(29) 



(30) 



/" rp( x >°) 2D 
Fimt(p(x,0)) = dx dz (p{x,0) - z) — - 

J Jr(x) V\ 



Remark : From the expressions (1211291 1 and ( [T9l , 

2 D(z) 

r(x) °{ z ) 

one can see that the case where D(p) = 1 and a(p) is a quadratic function of p, 

D(p)=l ; a(p) = 2Ap(B-p), (31) 
can be easily related to the SSEP, for which (see ifioll ) 

D(p)=l ; a(p) = 2p(l - p) . (32) 

In fact, if one makes the change of variable 

Bp ; j — > Bj 

one gets for the choice d3~Tl i that, both in the annealed and in the quenched case, 

p(\, p a , pb) = -^pssep (ABX, ^ ^) ■ (33) 
In the annealed case, where the exact expression of the SSEP is available (111213b . one gets 



Aiannealcd(A,/9 a ,Pb) = / .//.1<>U 



1+uie k 



where 



(34) 



Pa(e ABX -l) ^ Pb (e- ABX - 1) p aPb (e ABX - l)(e~ ABX - 1) 
« = B + B + W ' (35) 

In the limit A = B^ 1 — > 0, this gives /^annealed when a = 2p, i.e. in the case of non-interacting particles d49l that we will 
discuss in section[IV] Assuming that ( 1341351 ) remain valid for non-physical values of A and B, one would get /z annea i c d, without 
any further calculation, for the Kipnis Marchioro Presutti model flq, Il7l [l8l [l9ll where a = 4p 2 (in the limit B — > 0, A — > —2). 
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III. THE TIME REVERSAL SYMMETRY 

In this section we are going to see that the symmetry (0]i can be extended to more general diffusive systems. To do so, let us 
consider the difference T-^t, {p(x, 1)) — J-i n it.(p( x i 0))- Using ( fT8] l, one has 

/oc 
f(p(x, 1)) - f(p(x, 0)) - (p(x, 1) - p(x, 0))/'(r(x)) dx , 
-OO 

which can be rewritten as 

KnMx, 1)) - T iuit .(p(x, 0)) = / dr ^ dxf'(p(x, r)) 9p(x ' T) - H (p(x, 1) - pfo 0))f'(r(x))dx . 

J0 J -oc " T J -oc 

Then, using ( fT3] l, an integration by parts, and (fTTJ, one gets 

F init .(p(x,l)) - T init .(p(x,Q)) = j dr j°° dxf'(p(x, T ))^^-j(x, r) - j°° (p(x, 1) - p(x, 0))f'(r(x))dx 

'dr r dx 2D{P } X '\l )dpi ^ T) j(x,r)- r (p(x,l)-p(x,0))f'(r(x))dx. 



— oc 



cj(p(x,t)) dx 



This allows one to rewrite the last term in (f2Tb as 



^b- + ^P] 2 ^ i n it .(^l))-^ 1 ,(p(^0)) + r ip{xA) _ p{x , 0)) fmi dx+ ff dT dx^ + ^ 2 



2a 2 J_^ r " ' ^ " 2 77 2a 

(36) 

and therefore (|2"TT i becomes 



Aiannnealed(A) = max i h/ dx \fi{X, 1) - p{X, 0)\ 

P(x,T),j(x,T) 



2 

\ r , 3(x,T)l + (D(p(x : T))d xP (x,T)f \ 

J-oc 2a(p(x,r)) ' 



For the step initial density profile (1201 , one has (fTTT i 



f'(r(x)) = f'(p a )-6(x) I 2 -^-dp. 



in, 



°{p) 



One can then see in (l37l i that the initial time and the final time play symmetric roles : if one replaces {p(x, t), j(x, r)} by 
{p(a;, 1 — t), —j(x, 1 — r)}, (|37T > is left unchanged provided that A — > —A — ^^y-dp. (one has to use that, from the 

conservation of the total number of particles, [p(x, 1) — p(x, 0)]dx = J . Therefore /i a nncaicd satisfies 



Mannoalcd(A) — ^annealed f — A — / —^-dp) . (38) 

This is a generalization of (for the SSEP ( |32i >, 7J (p) = 1 and aUA = 2p(l — p), and ( |38l reduces to (01) and therefore 
shows that a version i32ll of the fluctuation theorem J40l l4ll l42l l43l H, l45ll holds for general diffusive systems with the step 
initial condition considered here. Although this initial condition is neither an equilibrium state, nor a non-equilibrium steady 
state, the time reversal symmetry d38l holds. We think that this is because, in the annealed case, the initial condition is in local 
equilibrium. 

One can repeat the same transformations in the quenched case. Due to the absence of J r i n it.(p(^, 0)) in (1291 ). one ends up 
with an expression where p(x, 1) and p(x, 0) do not play symmetric roles, so that /i quC nchcd does not seem to satisfy any kind 
of time reversal symmetry. 



Remark : By the same reasoning, one can show that the symmetry ( 1381 ) holds for initial conditions more general than 
the step initial profile. One can consider at t = an initial density profile 

r{x) = (1 - v(x))p a + v(x)pb , (39) 
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where v(x) is no longer a step as in ( f20b but could be a more general sigmoid function with v(— oo) = and v(oo) — 1. One 
can also replace the measure of the integrated current ( |T6b at the origin by its weighted average over space in a region around 
the origin : 

/oo 
dx w(x) [p(x, 1) - 0)] , 
-oo 

where w(x) is another sigmoid function. Then following exactly the same steps as in the derivation of (T37l > one gets 



Mannnealed(A) = max 

p(x,r),j(x,r) 



dx [p(x, 1) - p(x,0)} 



— oo 

1 

dr 



Xw(x) 



Pa 



dx 



j(x,r) 2 + (D( P (x,r))d x p(x,r)) 2 



dp 



2a(p{x,T)) 

from which one can see that the time reversal symmetry d38l remains valid if v(x) and w(x) are related by 



w(x) 



r(x) 



Dip) 
<P) 



dp 



I 



pb 



Dip) 
<r(p) 



dp 



(40) 



Remark : No time reversal symmetry seems to hold in the quenched case. However, if an additional symmetry (the particle- 
hole symmetry) holds, one can relate /x annca i C( j and p qucnc hcd- In Appendix iBl we show that, if Dip) and a{p) satisfy 



Dip) = £>(2p-p) 
aip) = ai2p-p) ' 

then the optimal profile p^ (aj, r) for the annealed variational problem (l2TT i when p a = P b = p is such that 

p^ix,r) = 2p-pM(x,l -t) . 

This implies in particular that p^ix,r = 1/2) = p and allows one to relate the optimal annealed (f2T~b and quenched 
profiles (see AppendixlBt. leading to 

Mqucnchod(A, p a = Pb = P) — /^annealed (A, p a = Pb = p) ■ 

For the SSEP d32l the particle-hole symmetry (|4TT > is satisfied, and therefore d43l > holds, for p = 1/2. Thus 
Mqucnchod(A, p a = Pb = 1/2) can be deduced from the exact expression J2I3I ). 



(41) 



(42) 



(43) 



IV. THE NON INTERACTING WALKERS 

The problem with the expressions (l2Tb or d29b is that it is very hard to solve the equations satisfied by the time dependent 
density and current profiles for general Dip) and aip). In this section, we solve the easy case of non-interacting random walkers. 

Let us consider non-interacting particles on an infinite one dimensional lattice. Each paticle on this lattice jumps at rate 1 to 
each of its neighboring sites, irrespective of the positions of the other particles. One can show (see appendixlct that in this case 

D(p) = l ; aip) = 2p ; f(p) = plogp- p. (44) 

Then d23l becomes 

dpix,r) d 2 pix,r) d 

dr dx 2 dx 

and the evolution equation (l24"i i of H becomes autonomous : 

8Hjx,T) _ d 2 Hjx,r) f dH{x,T) V 
dr dx 2 \ dx J 

It is easy to check that the general solution of ( 145|46t can be written as 

Hix, r) = K + lnG(a;, r) and r) = G(a;, r)Rix, r) , 



2 P (x,t) 



dHjx,r) 
dx 



(45) 
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where G antidiffuses and R diffuses : 

8G(x,t) 8 2 G(x,t) 8R(x,t) _ d 2 R(x,T) 

dr dx 2 dr 8x 2 

As the boundary condition (l25T l holds both for the annealed and the quenched case, one gets 

e A + 1 e A — 1 
G{x,r) = e -4± + e —±E 



E(z) = —= \ e~ u du. 

\/7T 



where E(z) is the error function 



In the annealed case, the boundary condition d26t becomes 

p(x,0) = r{x)e- X6{x) G(x,0) 
and, using (l44l >. the solution of (l45T l for this boundary condition is 



(47) 



p(x,r) 

Using ( 123124144b . one can show that 



PbS A + Pa , Pbe A - p a ( 



8H\ 2 8{Hp) 8 ( TT 8p 8H 
dx 



8t dx V"<9x ^ dx ^P q x 



H- 



G(x,t). 



dH 



Using this identity in d27l > and the fact that p and H have limiting values at ±oo, one gets 

/■oo />oo 

j«annnealed(A) = -Fia&.(p{ x i °)) + A / dx[p(x, 1) - 0)] - / dx[H(x, l)p(x, 1) - H(x, 0)p(x, 0)] 

JO J -oo 

Using ( 111118125b and fl26), one then has 



(48) 



Mannncalcd (A) = 

and, as f(p) = p log p — p, one gets 



dx[p(x, 0)f'( P (x, 0)) - f(p(x, 0)) - r(x)f'(r(x)) + f(r(x))} , 



Mannncalcd (A) 

This leads to 



dx[/3(a;,0) - r(x)] = p a - 



/^annncalcd (A) 



1 f° 
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-(f)] 


J —oo 






Pa(e x - 


1)- 







dx + pb- 



e~ x - 1 



1 -£ 



(f) 



(49) 



One can notice that this is just the limit of ( |2|3t when p a and p\> are small (at low density the exclusion rule in the SSEP can be 
neglected). One can also see by expanding (l48b in powers of A that in the long time limit 



(Q) Pa - Pb . (Q 2 )c ; Pa + Pb 
t Jtt t Jtt 



(50) 



In the quenched case, the boundary condition is 



p(x,r) = 



Pb 



instead of 

Pa , Pb~ Pa 



E 



Therefore the profile becomes 

G(x,t) 



2^ 



G(x,0) 
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Then, following the same steps as in the derivation of d48l l. one gets 



Mquonchcd(A) = Xj dx[p(x, 1) - p(x, 0)] - / dx[H(x, l)p(x, 1) - H(x, 0)p{x, 0)] 



which leads to 



/0 poo 
dxlogG(x,0) + pb dx\og[e~ x G(x, 0)] 
-oo JO 

Therefore 

Mquonchcd(A) = p a I dx log 
J —oo 

The expansion in powers of A leads to 



e x + 1 e A - 1 fx 
E 



2 V2 



Pb / dxlog 
h 



l + e~ A 1 - e~ A fx 
= + = E' 



(Q) Pa- Pb . (Q 2 )c . Pa+Pb 
t ~* 



which shows that the annealed ( T50b and quenched (l52l cases start to differ at the level of the variance of Q t . 



(51) 



(52) 



Remark : Taking the A — * oo limit of the expression (Bil l of ^ qU enched( A), one obtains 

Mqucnchcd(A) ~ TTPaA 3 ^ 2 , 

A — >oo o 

(with a similar result with pb replaced by p a and A by |A| for A — > — oo). Then, we can perform a Legendre transform to obtain 
the decay of the distribution of the integrated current Q t , as defined in (|7]), which yields 



Pro 



'Qt 




[ Vtq 3 ] 




x exp 

Q — »OG 


[ 12pgJ 



(53) 



This non-Gaussian decay is very reminiscent of the SSEP (0. In section [VTJ we will show that this type of decay is rather 
generic. 

Expression ( 1531 can alternatively be understood from (IC21 l. as the tail is dominated by the contribution of the first Q t particles 
at the left of the origin, that is : 



Pro 



Qt 



Qt 

J|exp 

i=l 



exp 



where we have used that the average distance between consecutive particles is l/p a . In the annealed case where the initial profile 
can fluctuate, the decay is slower, because the events which dominate have an initial profile where the Q t particles are arbitrarily 
close to the origin. 



V. ROTATIONAL SYMMETRY FOR THE SSEP 

In this section, we consider the MFT of the symmetric simple exclusion process (SSEP), for which d32l D(p) = 1 and 
(j(p) = 2p(l — p). The MFT then exhibits a remarkable symmetry : in the annealed case, this symmetry allows us to relate the 
generating functions of the integrated current Q t for different values of the initial densities p a and pb- This relationship takes the 
form of the single-parameter dependence (01 

Mannealed(A, p a , Pb) = F(u(X, pa, pb)) , 

with u> given by ((3). This to dependence was already derived by considering the microscopic dynamics of the SSEP in l29ll . 
Here, it is recovered by showing that, when uj(X, p a , Pb) — w(A', p' a , p' b ), an explicit transform relates the variational problems 
( t2Tb with parameters (A, p a ,Pb) and (A', p' a ,p' b )- 
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This transform is inspired by a known representation of the microscopic exclusion process |H1] in terms of spins : here, the 
equivalent of a global rotation of these spins will allow us to go from (A, p a , pb) to (A', p' a , p' b ). When p is expressed as an 
optimum (f22b over the two independent variables p(x, t) and H (x, i), one can introduce a "spin" variable, 



S+ =pe~ H 
S- =(1-/ 
S z = p- \ 

The bulk term in the variational problem ( 1221 can then be rewritten as 



S- =(l-p)e H , and the quadratic form S- S' = -(S+S'_ + SS' + ) + S% 



// 



Hd T p - d x Hd xP + p(l - P ){d x H) 2 = / / -Hd T p - d x S ■ dj . (54) 



The last term of this "action", —d x S ■ d x S, is clearly invariant under orthogonal transforms of S. Thus, starting from the 
optimal profiles (p, H) for a given set of parameters (A, p a , Pb), one can deduce sets of profiles (//, H'), obtained by performing 
an orthogonal transform on S, which satisfy the same bulk minimization equations (12311241 as (p, H). 

Therefore, for (//, H') to be the optimal profiles for other values of the parameters (A', p' a , p' b ), it is sufficient that they satisfy 
the corresponding boundary conditions : (125126b in the annealed case, and ( 125128b in the quenched case. 

Let us first look at these boundary conditions at r = 0, 1 for x — ► ±oo : 

p(— oo.t) = p a „ f p(oo,t) = Pb 

1 K ' ; p as well as > yy ' ; ' 



#(-oo,t) = \H(oo,t) = X ' 

which correspond to 5-00 = (p a , 1 — p a , p a — 1/2) and S+oc = (pbe~ x , (1 — pb)e x , pb — 1/2). Under an orthogonal transform 
on S, the scalar product of these vectors is necessarily conserved : 



i (p a (1 - p b ) e x + p b e- x (1 - p a )) + (^p a - ( Pb - = I + i 



with oj as defined in (O. Hence u = w' is a necessary condition for (p' , H') to be optimal for the set of parameters (A', p' a , p' b )- 
In order to explicitly check that one can indeed relate the optimal profiles when u> = u>', and to compare the corresponding 
generating functions, we will now express the optimal profiles (p, H) for (A, p a , Pb) in terms of the "reference profiles" (p, H) 
obtained for the SSEP at uniform density 1/2 : p a = f>b = 1/2. When uj — u, we reparametrize p a and pb in terms of two 
variables u and v : 

cosh it— 1 e~ v cosh u — 1 

Pa = t : and p b = , 

e A — 1 e A — 1 

so that 

lu = sinh 2 u and A = 2u . 
One can then check (after some algebra) that the mapping (p, H) — > (p, H), as 



P = ■ , a e H - u sinh _ S i n h A^ii ~ pe u-H sinh u±v _ ^ _ -x sinh u=v 

„i? = 1 , e"(e A -l)(e a -l) 

e H (e"-e")+e«(e" + «-l) 



(55) 



gives a solution of the bulk equations (I23l24l >. 

From the expression of e H , one can easily see that the final time boundary condition (l25T l, which is common to the annealed 
and quenched cases, carries over from H to H : 

= 2u0(sg) H(x,l) = X9(x) . (56) 

However, the initial-time boundary condition behaves differently in the annealed and in the quenched cases. In the quenched 
case, one would need that p(x, 0) = r(x) when p(x, 0) = 1/2 : this requires d55l l that H(x, 0) = X8(x), which is not expected 
to be satisfied as H (x, 0) is free under the quenched boundary conditions. Hence the condition d28| i does not carry over from 
(p, H) to (p, H), and d55l l does not lead to the correct optimal profiles in the quenched case. 
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On the other hand, the initial-time condition in the annealed case (l26l l is H(x, 0) = 2u6(x)+f'(p(x, 0)) — /'(1/2). Integrating 
the Einstein relationship (fTTT i for £) = 1, cr = 2p(l — p) leads to 



/'(r) = log and f(r) = r log r + (1 — r) log(l — r) . 

1 — r 



One can then check that ( |55t yields 



#(a;,0) = 2u0(z) +log 



p(s,0) 



H(x, 0) = A0(a;) + log - ^ °L - log ■ ^ 



l-p(x,0) °l-r(i) 



(57) 



(58) 



Therefore d55l l maps the optimal profiles for (A, 1/2,1 /2) to those for (A, p a , pt,) in the annealed case. 

This in turn allows us to relate the generating functionals /x annea i e d(A, 1/2, 1/2) and /i an neaied(A, p a , Pb) ■ taking into account 
the invariance of the bulk term, we obtain from (1221541 ) 



Pan. (A, Pa, Pb) - Pan. (A, 1/2, 1/2) 



dx [X(p(x, 1) - p(x, 0)) - X(p{x, 1) - p(x, 0)) 

Finit, (p(%, 0)) + f init. (p(x, 0)) 

dTdx[Hd T p — Hd T p] . 



Integrating by parts the last term and using j56ll58] ), this can be simplified to 

/w(A,p a ,p 6 )-M M .(A,l/2,l/2) = J ^ lo g /j (^o) V-l/2 0) + // dTdx iP d r H - pdrH] 

From ( f55l l. one can express pd T H — pd T H as a total derivative in terms of H : 

pd T H - pd T H = - JL fog [( e « _ e ")e A + e u (e u+1 ' - 1) . 

or L J 

Then, using the boundary conditions (l56l l and ( 158b as well as d55l l. we can evaluate d59l ) : we obtain 

1 - r(x) 1 - p(x, 0) _ (e u - e v )e S ( x <V + e u (e u+v - 1) 
l-p{x,0) 1-1/2 ~~ ( e « _ e^gif^o) + e «( e n+« _ 1) 

at each x, so that p an .(X, p a ,Pb) = Man. (A, 1/2, 1/2). 



(59) 



VI. BOUNDS ON THE DECAY OF THE CURRENT DISTRIBUTION 



In this section, we attempt to generalize the non-Gaussian decay (15153b of the distribution of the integrated current Q t during 
time t, 



Pro 



Q_t 



"Vtq 3 



to other diffusive systems. We had a = for the SSEP in the annealed case ll29ll . and a = for non-interacting particles in 
the quenched case d53l l. 

Here, we show that this form of decay holds, both in the annealed and quenched averages, when the following conditions are 
satisfied : 



D{p) = 1, 

f(p) < p + cfoiO < p < R, with a(p) = otherwise. 



More precisely, we set out to show that, when t — ► oo then 

.3 ! 



2PaV{Pa) Vt 



< — = log Pro 



< 



12(i? + c) 



(60) 



(61) 
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Let 



g(q) = lim —= log Pro 



In the MFT, g(q) is expressed as the optimum of a variational problem, like the current generating function //(A) (see ( 121129b ) : 

5(9)= , max {-^14.(^(^,0))- // drdx ^* 9 *^ X , (62) 

p(x,T),j(x,r) I J J 2a(p) J 

where the density profile p(x, t) is such that dx [p{x, 1) — p(x, 0)] = q, and where the current profile satisfies the con- 
servation law d x j + dtp = 0. In addition p(x,0) is free in the annealed case while it is constrained to be equal to 
r(x) = p a + (pb — p a )0(x) in the quenched case : hence 

^/quenched \0) ^ ^annealed 

Let us first obtain the lower bound in doTT i. Because of the variational formulation ( f62b (in the quenched case, Fi n it.(p( x > 0)) = 
0), one can bound g(q) from below by considering a particular profile (p(x, t),j(x, t)) leading to a total flux q. Here, we choose 
to move the segment [—q/p a , 0], which contains q particles at time 0, at constant speed v = q/p a from time to time 1, so that 
the total flux through during this time will be exactly q : this corresponds to 



q for - q(l - r)/p a <x< qr/p a : 

—d x p otherwise. 



Since p(x, r) = p a for —q(l — r)/p a < x < qr/p a , this leads to 



g(q) >~ I dr [ dx- 



-q(l-T)/p a 2 v(Pa) 2p a o( Pa ) ' 

which is the lower bound in doTb both in the annealed and in the quenched cases. 

The upper bound is obtained by noticing that, if a(p) = outside of [0, R] as in d60b . the fluctuation-dissipation relationship 
(1 1 01 lib implies that Finit.(p( x i 0)) diverges for p(x, 0) ^ [0, R]. From d60b , i.e. <j(p) < p + c, we then obtain 

9 (q)< max - [ [ drdJ j + ^ 



P (x,T),j(x.j) J J 2(P + C) 

where p(x, r) is such that < p(x, 0) < R and J °° dx [p(x, 1) — p(x, 0)] = q. The right-hand side of (?? ) is the maximum 
over p(x,0) of the g qU enched(9) f° r non-interacting walkers with initial density p(x, 0) + c : it is maximal, for q > 0, when 
p(x, 0) is equal to R for x > and for x < 0. This corresponds to the quenched, non-interactive case 031 at densities R + c 
and c, so that 

q 3 

9(q) < - J 



12(i? + c) 2 
which is the upper bound of ( |6Tb . 



VII. CONCLUSION 

In the present work, we have shown ( |23ll30l l how to implement the macroscopic fluctuation theory to study the fluctuations of 
the current of diffusive systems with a step initial density profile. We have argued that, depending on whether the initial profile 
can fluctuate or not, one has to perform an annealed ( 12 1 127b or a quenched average d28H30b . Using the structure of the equations 
to be solved in the MFT, we could obtain a simple relation ( 13 1133b between the generating functions of the current of the SSEP 
and of other models with a quadratic a(p) such as the Kipnis-Marchioro-Presutti model. Thus our solution l29ll for the SSEP 
determines the generating functions of the current for all these other models. We established in section [Til] that a time reversal 
symmetry ( 138139140b . which is a version of the fluctuation theorem for a non-steady state initial condition, holds in the annealed 
case. In section|lV]and in Appendix|C]we showed that the case of non-interacting particles can be solved both by a macroscopic 
and a microscopic approach. In section[V]we have seen that the u> dependence of the SSEP could be understood as a rotation 
invariance of the MFT and we have exhibited d55b how the optimal profiles are changed under these rotations. Lastly, in section 
I VII we have shown that the non-Gaussian decay (0 of SSEP is generic under some simple conditions on cr(p). 

The main difficulty that we could not overcome was to solve the equations ( T23ll26l28"l ) satisfied by the optimal p(x, r) and 
H(x, t), even in the case of the SSEP where the generating function is known. Even for large A, we were unable to solve them, 
which is why we could only get bounds on the decay of the distribution of the integrated current Q t in section[VT] Solving these 
equations, even in the large A limit, remains an open question. 
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APPENDIX A 

In this appendix, we first show, as in [8|], how the variational form $2l[ where one has to optimize over density and current 
profiles which satisfy the constraint (TT3T > can be replaced, using the Martin-Siggia-Rose formalism, by the expression (l22l where 
the profiles p(x, r) and H (x, t) do not satisfy any constraint. We then show that the optimal p(x, t) and H(x, r) are solutions 
of ( 1231241 ) with the boundary conditions ( 125126b . 

Let Pro(po(a;) — * pi(x)) be the probability of observing the rescaled density profile p\(x) at time t, starting from an initial 
profile po(x). Formally, it can be written (TT~5b as a functional integral over all the density and current profiles (p(x, r),j(x, t)) 
statisfying p(x,0) — po{x) andp(x, 1) = pi(x) : 



Pro(p (x) -4 pi(x)) ^ / VpVj 



J]_5(d T p + d x j 



exp 



V~t dxdr 



(j + D{p)d xP f 
Mp) 



where the constraint dot appears as a 6 function at each point (x, r). One can then use an integral representation for each of 
these S functions by introducing a new field H (x, r) : 

(j+D( P )d xP r 



Pro{p {x) -»■ Pl (x)) x / VpVjVHexp 



-Vi / / dxdr H(d x j + d T p) 



2a(p) 



One can integrate by parts J dxHd x j (this entails no boundary term as j is expected to vanish at ±oo) to express the right-hand 
side as 



VpVjVH exp 



v7 // ,/„/, (Hd T p + D(p) dxP 8 x H-°M(8 x H? + U + D(p)d^-.{p)d x H)- 



(Al) 



After a Gaussian integration over the currents j(x, t) we obtain Pro(po(x) — > pi{x)) as an integral over the two unconstrained 
fields p and H : 



Pro(po(x) pi(x)) x J VpVH exp - JJ dxdr \ Hd T p + D(p)d x pd x H 
Taking ( IA2b together with ( [ToT l and ( fT7b . one gets /Lt ann eaied(A) as a extremum over p and H : 



-(d x Hy 



A < annnealod(A) = max 
p,H 

which is d22l) . 



-FinMx,0)) + \ J dx \p(x,l) -p(x,0)] - JJ drdx [Hd T p + D{p)d xP d x H 



(A2) 



(d x H) 



One can then determine the equations satisfied by the optimal profiles for p and H by looking at the effect of a small 
variation, p(x, r) — > r) + r) and H(x, r) — * i?(a:, r) + SH(x, r) : after a few integrations by parts, one obtains 







dx8p{x, 0) 



5Fu 



X9(x)+H(x,0) 



dxSp(x, 1) [\6{x) - 1)] 



d T H + D( P )d 2 x H+^(d x H) 2 



Sp(x,0) 

+ JJ drdx SH(x, t) [-dr P + d x {D{p)d xP - a{p)d x H)] 

+ JJ drdx 6 P (x,r 

This yields the two bulk equations ( 1231241 satisfied by p and H at the optimum : 

( d T p = d x (D(p)d x p - <r(p)d x H) 

{ d T H = -D(p)d 2 x H - ^(d x Hf ■ 

The first of these equations is just the conservation law, d x j + d x p = 0, since, from ( lAlb . we have 

j = -D(p)d xP + a{p)d x H 



at the optimum. Using ( TT~8T > to express 



6p(x,0) ' 



we also obtain from (IA3b the boundary relationships 



H(x,l) = X9(x) 

H(x, 0) = X0(x) + f(p(x, 0)) - f(r(x)) ' 



which reduce to ( 1251261 1 by using dTTT >- 



(A3) 
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APPENDIX B 

In this appendix, we show that when D(p) and a(p) satisfy the particle-hole symmetry d4TT >. the optimal profile (assuming 
that it is unique) in ( f2TT > verifies ( |42T > when p a = Pb = P- This will allow us to relate the optimal profiles in the annealed and in 
the quenched cases and to obtain (1431 . 

First, when p a = pb(= p), the term proportional to f'{r(x)) in (|37| | vanishes due to the conservation of the total number of 
particles, so that fl37| > becomes 

Mannnoalod(A) = max <^ VI dx \fi(x, 1) - p{X, 0)J \0(x) 

p(x,T),j(x,T) I Z J-OO 

;r fe i(v)' + w,rM,r)n (B1) 



Mp(*.t-)) 

Moreover (|4TT> implies (see([T9b) that 

Therefore, if (p^ (x, r), j( a ) (x, r)) is optimal in dBll ). then (2p — p(°)(a;,l — t), j( a ) (x, 1 — r)) is also optimal and, if this 
optimum is unique, one gets d42i > 

p (a) {x,T) = 2p-p (a) (x,l -r) . (B2) 

Due to this symmetry, one can rewrite ( IB II ) as 

Mannnealed(A) = 2 max ( ^j^-£-El ^ _|_ f dx \f(x,l) - p]X6(x) 
p (x,t)J(x,t) I Z J_ oc 

1 „ f dx i(^)!±M4^1M^| . (B3) 

1/2 J-oo 2cr(p(x,T)) J 

with pto (a;, r = 1 /2) = p from dB2l . 

For the quenched problem, using the identity ( |36l l, the fact that the term proportional to f'(r(x)) vanishes, and that 
Fmit.(p( x > 0) = P) = 0' one can rewrite d29l > as 

{^"' (p(x 1)) /*oo 
; h / dx [p(x, 1) - p(x, 0)]A6»(x) 

1 dr f dx ^^(D^r^rm 
J-oo 2a(p{x,T)) J 

with the initial-time condition p(x, 0) = p. 

We see that (IB3t and ( IB4l i are identical except for the range of variation of r. This allows us to relate the optimal profiles in 
the annealed and the quenched cases by 



from which ( 1431 follows easily. 



APPENDIX C 

In this appendix, we first show why, for non interacting walkers on a one dimensional lattice as in section|IV] D(p), a(p) and 
f{p) are given by (l44l >. We then explain how d49b and ( Bit can be recovered by a microscopic calculation. 

Consider first a Id lattice of length L : a new particle is injected at rate a on site 1 and at rate 6 on site L. Each particle on 
site 1 is removed at rate 7 and on site L at rate i5. As the particles do not interact, the probability that a particle Tj on site i will 
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have escaped, after time t, into the right reservoir evolves according to 

T 2 -(l+ 7 )T i; 



~~dt 

dt 
dT L 

dt 



whose solution in the long time limit is 



T i+ i + Ti-i - 2T j; for 2 < i < L - 1 ; 
p + T L _ X - (1 + 0)T L , 



T,= 



It is easy to see that the contribution to Qt of the particles entering the system during the first time interval dt is 

(e XQt + dt ) = (e XQt ) (l + aTi(e A - 1)* + 5(1 - T L )(e~ A - 1)*) . 

Therefore 

1, , ^(e A -l) + |( e - A -l) 



which becomes for large L 



lim -log(e AQt ) = - . 



lim i]og(^> = ^(« A -D+^- A -l) 



with p a — — and p — j^. The expansion in powers of A (see (18191 leads to D(p) = 1 and a(p) = 2p, as in (PBi i. For these non 
interacting particles, the partition function Z(N, L) = L N /N\ so that 

f(p) = Plogp- P 

as in d44l > , and (fTTl i is verified. 

One can also see that at equilibrium, at density p, there is an invariant measure (the equilibrium) where the occupation numbers 
rii of the sites are independent random variables distributed according to a Poisson distribution 



Pro(n) 



Z(N -n,L- l)Z(n, 1) _ p n e~ p 
Z(N, L) ~ n\ 



(CI) 



Let us now consider non-interacting particles on an infinite one dimensional lattice. Each particle jumps at rate 1 to each of 
its neighboring sites. The probability P, j (t) that a particle initially at position i will travel a distance j — i is given, for large t, 
by 



1 U-i 

-e 44 



/Ant 



(C2) 



The contribution of a particle initially located at site i to e y * is 

*< = 1 + (1 - 0O(e A - 1) + ^( e_A 

where (?j = 1 if i > 1 and 9i = if i < 0. In the long time limit, this becomes 

1 + £(_*) 

$i = 1 + (1 - 0;)(e A - 1) + 6 l (e- x 



j<0 



1)- 



where £7 is the error function defined in ( |47T >. Therefore, for a given initial condition where the occupation numbers n, of all the 
sites are specified, one gets 



istory 



exp 



^ rii log $i 
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The rii are distributed according to a Poisson distribution (ICU of density p a on the negative axis and pi, on the positive axis. 
Averaging over the rii (i.e. over the initial conditions) leads to ( |49] l in the annealed case and to d5TT i in the quenched case. 
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